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We consider the free energy W[J] = W^(H) of QCD coupled to an external source J^(x) = 
Hn cos(k ■ x), where ff£ is, by analogy with spin models, an external "magnetic" field with a 
color index that is modulated by a plane wave. We report an optimal bound on Wk(H) and an 
exact asymptotic expression for Wt(H) at large H. They imply confinement of color in the sense 
that the free energy per unit volume Wk(H) /V and the average magnetization m(k,H) = k dW £^ 
vanish in the limit of constant external field k — > 0. Recent lattice data indicate a gluon propagator 
D(k) which is non-zero, D(0) ^ 0, at k = 0. This would imply a non-analyticity in Wf,{H) at 
k = 0. We also give some general properties of the free energy W(J) for arbitrary J(x). Finally 
we present a model that is consistent with the new results and exhibits (non)-analytic behavior. 
Direct numerical tests of the bounds are proposed. 
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1. Introduction 

Recent numerical studies on large lattices of the gluon propagator D(k) in Landau gauge in 
3 and 4 Euclidean dimensions, reviewed recently in [R1], yield finite values for D(0) ^ |Q] - [f7|], 
in apparent disagreement with the theoretical expectation that D(0) = 0, originally obtained by 
Gribov [|p, and argued in [|9p. The argument [Q] which leads to D(0) = 0, relies on the hypothesis 
that the free energy W(J) in the presence of sources / is analytic in J at low momentum k. That 
hypothesis should perhaps be dropped in view of the apparent disagreement with the lattice data. 
This is of interest because a non-anal yticity in the free energy is characteristic of a change of phase. 

The free energy W(J) enters the picture because it is the generating functional of the connected 
gluon correlators. In particular the gluon propagator is a second derivative of W (J) at J = 0, 



at , ^ 8 2 W(J) 



(l.i) 

7=0 



8J«{x)8J b v (y) 

The free energy W(J) in the presence of sources J is given by 

expW(J) = (exp(/,A)) 

= / dA p(A)exp(/,A), (1.2) 

JQ. 

where p, v are Lorentz indices, and a,b are color indices, and 

(J,A) = J d D xJ b il {x)A b il {x). (1.3) 

The integral over A is effected in Landau gauge dpAp = 0, and the domain of integration is re- 
stricted to the Gribov region Q., a region in A-space where the Faddeev-Popov operator is non- 
negative, M(A) = —d jl D jl (A) > 0. We use continuum notation and results, but we have in mind the 
limit of lattice QCD in the scaling region, that is gauge-fixed to the Landau (or Coulomb) gauge by 
a numerical algorithm that minimizes the Hilbert norm squared [|A|| 2 , and thereby fixes the gauge 
to the interior of the Gribov region. The vector potential, given by A(x) = gA psrt (x), is unrenormal- 
ized, and has engineering dimension in mass units [A(x)] = 1 in all Euclidean dimension D, while 
[H] = D — 1 . (Our results also hold in the Coulomb gauge at fixed time, in which case D is the 
number of space dimensions.) The density p(A) is a positive, normalized probability distribution 
with support in the Gribov region £2. Because there are Gribov copies inside Q., p(A) is not unique 
and, in general, depends on the minimization algorithm. 
We consider a source that has the particular form 

/J(x) =#Jcos(fcci), (1.4) 

where we have aligned the 1-axis along k, so the free energy 

e\pW k (H) = ( exp[ J d D x #Jcos(fcci)Aj(x)] ), (1.5) 

depends only on the parameters k and H b . This is sufficient to generate the gluon propagator for 
momentum k, 

D°?(k)=2 d ' W f\, (1.6) 

1 J 
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where 

, x W k (H) 

w k {H) = -*±-L (1.7) 

is the free energy per unit Euclidean volume. Because A ll (x) is transverse, only the transverse 
part of H is operative, and we impose k^H^ = 0, which yields Hf = 0, and we write Hf, where 
i = 2,... D. By analogy with spin models, Hf may be interpreted as the strength of an external 
"magnetic" field, with a color index b, which is modulated by a plane wave cos(fcti). (This external 
magnetic field Hf, with color index b, should not be confused with the Yang-Mills color-magnetic 
field Ft) 

A rigorous bound for Wk(H) on a finite lattice was given in [^] which holds for any (numerical) 
gauge fixing with support inside the Gribov region O. One can easily show that in the limit of large 
lattice volume V, and in the continuum limit, this implies the Lorentz-invariant continuum bound 
in D Euclidean dimensions, 

w k {H) < (2Dlc 2 ) 1/2 \H\, (1.8) 



where \H\ 2 = L^,*(^) 2 - A model satisfying the bound (1.8) was recently exhibited in [ |10| ] 



More recently, a stricter bound for Wk{H) at finite H was obtained [|1T|], that also holds for any 
(numerical) gauge fixing with support inside the Gribov region O, 

w k (H) < 2~ x l 2 k it[{H a H a fl 2 ]. (1.9) 

Here H a H a is the matrix with elements HfHj. It has positive eigenvalues, and the positive square 



root is understood. This bound is stricter than the old bound (1.8). It is in fact optimal for a 



probability distribution p (A) of which it is known only that its support lies inside the Gribov region. 



Expression ( J1.9| ) also provides the asymptotic form of w k (H) at large H, and infinite Euclidean 
volume V JlT| ] for any numerical gauge fixing with probability density p(A) with support that 
reaches all boundary points of £2, (but which may vanish on the boundary, p(A) = for A € d£l) 

w k ^(H) = 2- l ' 2 k tr[(H a H a )^ 2 }. (1.10) 

Either bound yields in the zero-momentum limit 

w (H) = hmw k (H)=0. (1.11) 

k— >0 

As discussed in [Q], this states that the system does not respond to a constant external color- 
magnetic field no matter how strong. It is a consequence of the proximity of the Gribov horizon 
in infrared directions. We shall return in the concluding section to the physical implications of this 
result for confinement of color. 

If w k (H) were analytic in H in the limit k — > 0, eq. ( |1.11| ) would imply that all derivatives of 



the generating function Wq(H) vanish, including in particular the gluon propagator (L6) at k = 0, 
D(0) = 0. However, as noted above, this disagrees with recent lattice data which indicate a finite 
value, 0(0) ^ 0, in Euclidean dimensions 3 and 4. If this is true, then w k (H) must become non- 
analytic in H in the limit k — >■ 0. In order to get some insight about this, we examine the behavior 
of an improved model that has the exact asymptotic behavior ( |1.10| ). 
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2. General properties of free energy in QCD 

The proof of the above results relies on properties of W (J) that hold for arbitrary J" (x) that are 



notable for their generality and simplicity [T\\{. The asymptotic form of W(J) at large / is given by 



Was (J)=]im^l, (2.1) 

A— s-oo A, 

and because of the convexity of W(J) we have the bound 

W(J)<W U (J). (2.2) 

Moreover the asymptotic free energy is given by 

W n (J)=mstx Aeda (J,A), (2.3) 

where dQ. is the boundary of the Gribov region, known as the Gribov horizon. Let this boundary 
be described by the equation h(A) = 0, where h(A) is the so-called "horizon function". Then, by 
the Lagrange multiplier method, the asymptotic form of the free energy is given by 

W as (J) = (J A*) (2.4) 

where A* = A* (J) minimizes 

I(A) = (J,A)-Xh(A), (2.5) 
and A is a Lagrange multiplier. Thus A* (J) is the solution of 

and A is determined by 

h(A*) = 0. (2.7) 

3. Improved model 



The model is defined by the expression for the free energy 

k 2 H a H a \ 1/2 



^W=^){tr[(/ + - W ) -/J-trln^^^^) +2^/]}, (3.1) 

where > is an as yet undetermined function, and H a H a is the matrix with elements HfHj, for 
/, j = 2, ... D. This model possesses the following desirable features [pi]: (i) It satisfies Wjt mo d(0) = 
0, which is correct at H = for a normalized probability distribution J dA p(A) = 1. (ii) It has the 
asymptotic limit 

w ku (H) = lim W ^ mod(At//) = 2- l / 2 k tx[(H a H a W 2 }, (3.2) 

that is correct at large H for any numerical gauge fixing that is strictly positive in the interior of the 
Gribov region CI. (iii) It satisfies the optimal bound 



Wk. 



mod(H)<2- l / 2 kti[(H a H a ) 1 / 2 }, (3.3) 
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which implies that the generating function vanishes at k = 0, wo, mo d(#) = 0. (iv) The matrix of 
second derivatives is positive in the sense that 

d 2 , k d (H) 
1 dHfdH b J ~ 

1 J 

holds for all vf and Hf, as required for this matrix to be a covariance. 

Because of the property wo, m od(#) = 0, there must be some non-analyticity if, as indicated by 
numerical calculations, the gluon propagator D(k) at k = is positive D(0) > 0. It is instructive to 
see what kind of analyticity this would be in our model. Let X{H) > be the largest eigenvalue 



of the matrix HfHj. Inspection of (3T) shows that Wk, m od(H) is analytic in H inside a radius of 
convergence 

H H ) = ^-- (3-5) 



Moreover from (3. 1 ) we have at small H, 



w k , mod (H) = —HfHf + g(k)0[k 4 H 4 /g 4 (k)]. (3.6) 

For the gluon propagator D(k) ~ d g^ffi ~ k 2 /g(k) to be finite at k = 0, as suggested by the 

lattice data, we must have g(k) = const k 2 near k = 0. In this case the coefficient of the H 4 term is 
of order l/k 2 , which diverges as k — > 0, as do all higher order coefficients. Moreover the radius of 
convergence of the series expansion of w^mod^) is h(H) = 0(k 2 ), which vanishes like k 2 . 

Suppose that g(k) has a power law behavior g(k) ~ k v at k = 0. Then the radius of convergence 
behaves like X(H) ~ k 2v ~ 2 , which approaches with k for v > 1. The gluon propagator behaves 
like D{k) ~ k 2 ~ v , and Wk, m od{H) is non-analytic in H at k = when the propagator has a power 
law D(k) ~ k p with p < 1. Gribov's original calculation gave D(k) ~ k 2 /m 4 which corresponds 
to g(k) = 0(m 4 ), and Wk, m od{H) is analytic in H at k = 0, with a radius of convergence X(H) = 
0{k- 2 ) -> oo for k ->■ 0. 



4. Conclusion 

By analogy with spin models, we define, for each momentum k, the analog of the bulk mag- 
netization in the presence of the external "magnetic" field Hf, 

which describes the reaction of the spin system to the external color-magnetic field. Its physical 



meaning in gauge theory is apparent from (|L5j) which yields, 

M b ^(k,H) = (j d D x cos (feci )A b i {x)) H 

= (\/2){a b (k) + a h {-k)) H . (4.2) 
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Thus the "bulk magnetization" is in fact the k-th fourier component of the gauge field in the pres- 
ence of the external magnetic field. We also define the magnetization per unit (Euclidean) volume 



given by 



mmH) = *SMl, (4.3) 



«*.») = (4.4, 



The asymptotic free energy ( |1.10| ) determines the asymptotic magnetization per unit volume 
at large H, 

= 2- 1 / 2 k[(H b H b )- 1 / 2 ] ij Hj. (4.5) 
Its magnitude is given by (m" as ra- ! as )(&,//) = k 2 /2, and we obtain the simple formula 

lim (mfmf)(k,H) = k 2 /2, (4.6) 

ff-*x> 

which holds for any numerical gauge fixing with support extending up to the boundary of the 
Gribov region £2. 

We arrive at the remarkable conclusion that in the limit of constant external magnetic field, k — > 
0, the color magnetization per unit volume vanishes, no matter how strong the external magnetic 
field, 

lim ]immf(k,H)=0. (4.7) 

Thus the system does not respond to a constant external color-magnetic field. In this precise sense 
the color degree of freedom m b (k,H) = Jy(a b (k) +a h i {—k))H is absent at k = 0. This conclusion 
holds whether or not the free energy Wk(H) is analytic in H in the limit k — > 0. Lattice data would 
indicate that it is not analytic. Besides reporting this result, we have presented a model, defined in 
(3.1), which saturates the asymptotic limit ( 1.10| ), and exhibits confinement of color. As we have 



seen, W^ m od(H) may be either analytic in H, or not, at k = 0, depending on the behavior of g(k) 
at k = 0, but in either case, the conclusion stands, that the constant color degree of freedom of the 
gauge field is confined. 

Equations (Q and ( p0| ) may be checked numerically, at least in principle, by using the 
formula expWk(H) = (exp[/ d D x H b cos(kx\)A b (x)]) to make a numerical determination of the 
generating function itself. For large values of H this may fluctuate too wildly. Alternatively one 
may measure the magnetization from the formula M b (k,H) = (J d D x cos(fcci)Af (x))h, where the 
source term H b cos(kxi)A h (x) is included in the action that one simulates. This requires simulating 
the theory fixed in the Landau gauge instead of generating an ensemble from the gauge-invariant 
Wilson action then gauge fixing. It may be convenient to do this by numerical simulation of stochas- 



tic quantization Q13Q because that avoids calculating the Faddeev-Popov determinant explicitly. 

Finally we wish to emphasize the generality and simplicity of the results on W(J) for arbitrary 
J(x) that are presented in sect. 2. 
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